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Abstract 

We consider the problem of metastability in a probabilistic cellular automaton 
(PCA) with a parallel updating rule which is reversible with respect to a Gibbs measure. 
The dynamical rules contain two parameters (3 and h which resemble, but are not 
identical to, the inverse temperature and external magnetic field in a ferromagnetic 
Ising model; in particular, the phase diagram of the system has two stable phases 
when (3 is large enough and h is zero, and a unique phase when h is nonzero. When 
the system evolves, at small positive values of h, from an initial state with all spins 
down, the PCA dynamics give rise to a transition from a metastable to a stable phase 
when a droplet of the favored + phase inside the metastable — phase reaches a critical 
size. We give heuristic arguments to estimate the critical size in the limit of zero 
"temperature" {(3 — > oo), as well as estimates of the time required for the formation 
of such a droplet in a finite system. Monte Carlo simulations give results in good 
agreement with the theoretical predictions. 
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1 Introduction 



Metastable states are ubiquitous in systems undergoing first order phase transitions. During 
their lifetime (which can be very long indeed) these states are practically indistinguishable 
from equilibrium states. Nevertheless, they cannot be described in the framework of the 
equilibrium Gibbsian formalism [I], 0, |!|. Their analysis in terms of dynamical models has 
led to a deeper understanding of metastability by providing detailed descriptions of the 
"escape routes" from metastable to stable states in certain idealized limiting situations. 

Following earlier work on systems with long range interactions @, [1], Q, the pathwise 
approach to metastability was introduced in ||. It was then used in and || to study 
rigorously the escape from metastability in the Ising model with nearest neighbor interactions 
and a small external magnetic field, evolving via Metropolis Glauber dynamics in a finite 
periodic domain, in the limit of temperature going to zero. The effect of the boundary 
conditions on the exit path from the metastable phase was analyzed in [[J. Metastability for 
spin systems with different interactions was investigated in [|IU| and in f|TTJ the problem was 
considered in a more general context. 

The case of finite temperature, infinite volume and external magnetic field going to zero 
was studied in [|l!| ; this situation, very interesting from the physical point of view, is math- 
ematically much more complicated than the zero temperature limit. The finite temperature 
case also has been studied by means of Monte Carlo simulations (see e.g. and by 

transfer- matrix methods [jHfl ; a clear discussion of these results can be found in ||15|| . 

In all the above works (except for that of Penrose-Lebowitz dealing with deterministic 
continuum systems) the spin systems evolve according to a stochastic continuous time or 
serial dynamics, for which at most one spin of the system is updated at any time. In this 
paper we investigate metastable behavior in systems with parallel evolution, i.e., in which 
all the spins of the system are updated simultaneously, at integer times t = 1,2,3,.... 
In particular, we are interested in finding how the escape time and escape path from the 
metastable phase are influenced by the parallel dynamics. A natural setting for this question 
is that of probabilistic cellular automata (PCA), specifically those whose stationary measures 
are Gibbs states of Ising models with short range interactions. 

PCA were first studied in the Soviet literature of the early seventies [IB] and since then 
have been applied in many different contexts; in particular, their connections with statistical 
mechanics were investigated in [T7|, [13]. In this paper we will consider a PCA for which 
the dynamics depends on two parameters, f3 and h, and which has the property that its 
stationary states are Gibbs measures for a certain Hamiltonian H(/3,h). Here (3 plays the 
role of an inverse temperature and h that of an external magnetic field, but the coupling 
constants in H depend on (3 in a complicated way; in particular, H(j3,h) ^ (3H{h). As in 
the standard Ising model we have that when (3 is large enough and h is zero there exist 
two different stationary Gibbs measures for the PCA, characterized by nonzero average 
magnetizations ±m*, while for h ^ there is a unique stable phase. We then pose the usual 
question of metastability: if at large (3 the system is prepared in the minus (plus) phase and 
the magnetic field is chosen positive (negative) and small, how does the system reach the 
stable phase? For definiteness we will always consider the escape from an initial all minus 
phase. 
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In spin models with continuous time dynamics an important role in this transition is 
played by the stable configurations 0, f| f| [K|, which are fixed points of the evolution in 
the limit of zero temperature. For example, a rectangle of pluses of width greater than one 
inside a sea of minuses is a stable configuration for the nearest neighbor Ising model with a 
small positive external field. The tendency of such a rectangle to grow or to shrink by the 
repeated addition or loss of single sites, as a function of its size, yields the behavior of the 
exit from the metastable state. In the PCA discussed below all configurations are accessible 
at a single updating. Nevertheless, we will argue that, for large (3, the only configurations 
relevant for the description of the exit from the metastable phase are those in which the plus 
phase is inside well separated rectangles in the minus sea. We then discuss quantitatively the 
growth and shrinkage of such droplets and compare our theoretical prediction with results 
of Monte Carlo simulations. 

In Section [| we define our model and show that it undergoes a phase transition at 
low temperature, and in Section |3] describe the specific model on which we will focus in 
the balance of the paper. We discuss our heuristics on the critical behavior of droplets in 
Section f|, and compare theoretical and Monte Carlo results in Section |5|. Section |5] is devoted 
to some brief conclusions. 



2 Description of the General PCA Model 

Let A be a d- dimensional torus containing L d lattice sites, i.e., A C Z d is a cube containing 
L d points and having periodic boundary conditions. At each site igA there is spin variable 
o~(x) = ±1; the space {1, — 1} A of configurations is denoted by Q. 

To define the dynamics of the model we introduce the discrete time variable n = 0, 1, ... 
and denote by a n the system configuration at time n. All the spins are updated simultane- 
ously and independently at every unit time; the conditional probability that the spin at site 
x takes value r at time n, given the configuration at time n — 1, is 

Prob((x n (x) = r|a n _i) = p x (r\cr n ^i) = ~ 1 + rtanh/?(^ K(x - y)a n -i(y) + hj . (1) 

Thus the time evolution is a Markov chain on Q with nonzero transition probabilities Pa (77 1 a) 
given by 

Pk(vW)= I[pMx)W), v^ea (2) 

xeA 

The coupling is of finite range (K(z) = if \z\ > Zq, with z < L and typically zo <C L) and 
the coupling constants K{z) will be held fixed throughout our discussion. The parameters 
(3 and h play the role of inverse temperature and external magnetic field respectively, as 
discussed above. Note that for large \h\, cr n (x) = sgn(/t) with high probability, while for 
large positive /3, cr n (x) = sgn[X^ K(x — y)cr„_i(y) + h] with high probability. 

We say that a probability measure p(a) on the configuration space Q is stationary for the 
PCA if and only if it remains invariant under the dynamics, i.e., iff J2a Pa{v\< j )p{ <t ) = Pi 7 ])- 
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By the general theory of Markov processes there exists, for any (3, h, and A, a unique 
stationary measure z/f ' h for the PCA. We say that the PCA is reversible with respect to a 
measure p iff 

P A {v\a)p{a) = PMV)P{V) , Vex, 77 e CI (3) 

Summing (|3|) over a shows that any p satisfying (|3|) is stationary for the PCA; the converse 
need of course not be true. It is however straightforward to check that if K(z) = K(—z) 
then the process @ is reversible with respect to the measure 

u^ h (a) = Z- 1 J] e^W cosh[/3(]T K(x - y)a{y) + h)) , (4) 

xeA y 

where Z is a normalization constant. To see this one simply notes that p x (rj(x)\a) can be 
written as 

/ / \ 1 \ 1 exp[(3hrj(x) + (3J2 y K(x -y)7](x)a(y)} 

Px(v{x)W) = - r / . / \ I 7 1 • (5) 

2 cosh/3[Y iy K(x - y)a(y) + h\ 

The measure (|4]) must of course be the unique stationary measure referred to above. From 
(^) it is clear that v^ h (a) is a Gibbs measure for a Hamiltonian H(j3,h) with (generally 
many spin) interactions of finite range, which by our assumptions are independent of A: 



cosh/5 (j2K(x-y)a(y) + hj 



(6) 



Hence taking the limit A /* Z d yields a Gibbs measure v@' h for H which is stationary for 
the PCA on Z d , defined by the natural extension of the Markov process (Q) to Z d . 

The stationary measures for the infinite volume PCA need of course no longer be unique. 
It is known in general, however, that if one stationary translation invariant (or periodic) 
measure is Gibbsian then all such measures are Gibbsian for the same Hamiltonian ||T9fl . 
Hence to find all translation invariant stationary states of our PCA we need only investigate 
translation invariant Gibbs states for H(/3, h). Such an investigation begins with the ground 
states of the Hamiltonian. For the model considered here it is easy to see from (||) that if 
K(z) > for all z, and if the set of K{z) which are nonzero is not chosen in a very special 
way, then for h = there are exactly two ground states of H(/3, h), +1, in which a(x) = 1 
for all x, and —1, in which a(x) = —1 for all x, while for h 7^ there is only one ground 



state. It then follows from the Pirogov-Sinai theory [20] that for d > 2 and j3 sufficiently 
large there will be in general two extremal translation invariant Gibbs measures for h = 
and a unique such measure for h 7^ 0. By the argument above, the same conclusion holds 
for stationary states of the infinite volume PCA. We are thus exactly in the same set up as 
in the familiar ferromagnetic Ising model. We remark that although we are not dealing here 
with pair interactions, or even exclusively ferromagnetic interactions, it is easy to see that 
the measures (g) satisfy the both the FKG |2T| and GKS p2j inequalities. 

We may now pose the paradigm question of metastability: if we prepare the system in 
the starting configuration —1 and take h to be small and positive, how quickly and in what 
manner does the PCA reach its stationary measure? We want to answer this question in 
the limit of /3 — > 00, with A and h fixed, in which the stationary state is +JL, and hence 
may formulate the first part of the problem as that of estimating the first hitting time 
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r+i = inf{?T, > : a n = +1}, in the limit (3 — > oo, when the system is prepared in o~o = — 1. 
To answer the second part of the question, we need to describe the path that the system 
follows to reach +1; typically, such a path will involve the necessity of passing through one 
of a small number of critical configurations. 

As in the case of continuous (or serial) dynamics, the first problem is to understand 
the behavior of the "stable" configurations, that is to estimate the probability that a stable 
configuration grows or shrinks. Rather than discussing this problem in general we shall now 
focus on a specific model. 



3 A Special Model 

For the rest of this paper we will focus on one special model from among those specified by 
([!]): the two-dimensional model (d = 2) with K(z) = 1 for z G A and K(z) = otherwise, 
where A = {0, ±ei, ±e 2 } is the set consisting of the origin and its four nearest neighbors. 
Thus the probability distribution of the spin a n (x) is determined by the spins at time n — 1 
at the five sites in a cross centered at x. According to (|), the stationary measure of 
this system will then be 

4> h (a) = Z-\(3; A) exp[- £ U x (a; (3, h) + (3h"£ a x ) , (7) 

where U x (a;j3,h) = Uo(r^ x a; (3, h) with t x the shift operator (with periodic boundary con- 
ditions on A) and 

U (a;(3,h) = - £ J lA 0,h)a(A) = -logcosh(/? £ a{y)+(3h), (8) 

AcA yeA 

with u(A) = UytA&iy) for any A C A. The six coefficients J\a\(P, h) are determined by the 
six values which the J2 y eA a {y) can take. For h = only even values of \A\ occur, and we 
find 

U (a;(3,0) = -J (/3,0)-J 2 (/3,0) £ a x a y - J A ((3,0) £ a(x)a(y)a(z)a(w) , (9) 

{x,y}GA {x,y,z,w}cA 

with 

j (/5,0) = -3-log[(cosh5/3)(cosh3/3) 5 (cosh/3) 10 ] > , (10) 
16 

MP, 0) = h log[(cosh 5(3) (cosh 3/3)/(cosh (3) 2 ] > , (11) 
lo 

J 4 (/?,0) = -^-logKcoshS^^osh^V^oshS/^) 3 ] < 0. (12) 
16 

The pair interactions are thus ferromagnetic while the four-spin interactions are not, so the 
usual conditions for GKS inequalities are not satisfied. 
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4 The Time Evolution 



We shall first give a heuristic argument showing that the important configurations for exiting 
from the metastable state in the PCA are, as for the usual Glauber dynamics, isolated 
rectangles of pluses of minimum width two. We shall then describe, again on a heuristic 
level, the growth and shrinkage of one such rectangle. We find a critical value l*(h), for 
h < 1, for the length I of the smaller side of the rectangle such that, in the limit (3 — > oo, all 
rectangles with / < l*(h) will shrink to zero (except for some special values of h, for which 
the condition is I < l*(h) — 1) while those with I > l*(h) will grow, resulting in an escape 
from the metastable state. 

Let us begin by comparing Glauber dynamics — realized via the Metropolis algorithm, as 
is usual in questions of metastability — with the PCA dynamics considered above, focussing 
on differences which are relevant when h is small and (3 is very large. The former, for a spin 
system with Hamiltonian H = H(h) and inverse temperature f3, proceeds by spin flips at 
single sites, with the rate c(x; a) at site x in configuration a given by 

c(x-a) = l 1 > XH(a*)<H(a), f ) 

l ' 1 \expP[-H(a x )+H(a)}, if H(a x ) > H(a), 1 } 

where o~ x is the configuration obtained from a by flipping the spin at site x: 

**(v)={7l y) ' T y z x ' (14) 

l o"(2/), for y ^ x. y ' 

Since c(x; a) depends only on (3[H{a x ) — H(a)} and is independent of (3 if H(o~ x ) < H(a), 
the dynamical landscape is determined entirely by the function (/3H)(a), and the stable 
configurations, i.e., those invariant under the dynamics in the limit f3 — > oo, are the local 
minima of H. 

In contrast, the PCA dynamics permits transitions from one configuration to any other 
in a single updating; we will see, however, that this distinction will play only a minor role in 
the analysis of metastability. Recall that the probability of a transition from a to another 
configuration 77 in one time step is given by the product of the probabilities of spin flips 
at sites where a and n differ with the probabilities of non-flips at sites where they agree. 
Probabilities of all possible single site flips are shown in Figure [I]; it is clear that, at large 
(3, certain flips are almost sure to take place, while all others have exponentially small 
probability. Thus from an arbitrary initial condition we expect a very rapid evolution to a 
stable configuration, with further change taking place on an exponentially slow time scale, 
and involving primarily flips of single spins — in fact, on six different exponential time scales, 
well separated for h < 1 and very large (3, corresponding to the six slow spin flip processes 
of Figure |I[ The parallelism of the dynamics is of relevance during the first, rapid, phase, 
but the analysis of metastability involves the slow, essentially serial, second phase — although 
in certain cases we must consider the effect of a small number of unlikely single-site events 
occurring simultaneously. 

A second difference is that, although in the PCA dynamics transitions which lower the 
energy are generally favored over those which do not, the probability of a transition from a 
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to r] is not specified entirely by the energy difference H(rj; j3, h) — H(a; (3, h). In particular, 
for the specific model introduced in Section the single-step probability of flipping a spin 
which agrees with two of its nearest neighbors is exponentially small, even when such a flip 
is (energetically) favored by the magnetic field, and hence there are pairs of configurations a 
and 77 which differ at a single site but are such that the probability of jumping between them 
(in either direction) goes to zero as (3 — > 00. This is illustrated in Figure [I]. Consequently 
there are many more stable configurations for the PCA than for the Glauber dynamics. In 
fact it is easy to see that any configuration in which the value of the spin at every site agrees 
with that of at least two of its neighboring sites is stable. 

Despite the large number of stable configurations, however, it is rectangular droplets 
which are important for exit from the metastable state, due to the effect of the most rapid of 
the "slow" single flip processes of Figure |l|. We formalize this as follows. Let Q be the set of 
configurations in which every plus spin agrees with at least two of its nearest neighbors, and 
define T : Q — > Q so that Ta is the configuration obtained from a by flipping all the minus 
spins with at least two pluses among their nearest neighbors. For any a £ Q the sequence of 
configurations T k a, k = 0, 1, ... , is nondecreasing, in the sense that (T k+l a)(x) > (T k a)(x) 
for all x, and hence must reach a fixed point er*, in which the set of plus spins forms well 
separated rectangles (or bands around the torus) inside the sea of minuses. Moreover, if we 
take a as the initial condition cq of the PCA dynamics, and let E a be the event that for 
some n, a n = a* and the sequence (cr m ) is increasing for < m < n, then it is clear from 
Figure ^ that 

lim Prob£ CT = 1, (15) 

and that the time to reach a* is typically of order exp 2/3(1 — h). Thus the path for escape 
from metastability must pass through configurations in which all the pluses are inside well 
separated rectangles. 

We do not attempt to discuss the most general situation but instead consider the fate 
of a single rectangle; we expect that, as for Glauber dynamics 0, this is the key element 
in an analysis of metastability. Let us consider, then, a configuration r] for which all spins 
inside a rectangle of sides I and m are up and all other spins are down, and suppose for 
definiteness that / < m. We will say that such an I x m droplet is supercritical if, starting 
from rj, the system will reach the configuration +1 before it reaches —1, with probability 
which approaches 1 in the limit (3 — > 00; the droplet is subcritical if the reverse is true. 
We will argue heuristically that if I < 2/h then the droplet is subcritical and if I > 2/h 
the droplet is supercritical, while if 2/h is an integer and I = 2/h then the droplet may 
either shrink or grow. The critical length is the smallest integer such that the droplet is 
supercritical if I > l* h ; thus l* h = [2/h\ + 1. 

A very rough estimate of 1^ may be based on energy considerations, with the assumption 
that if the system starts from a rectangular droplet then the next droplet reached will be 
one with lower energy. For example, if rj contains an I x I droplet and e(Z) = H(rj), one may 
approximate e(l) at very large f3 by writing log cosh f3x ~ (3\x\ in (U). In this approximation, 
e(l) is a parabola whose maximum is achieved at I — 2/h, supporting the result l* h = [2/h\ +1 
described above. 

For a correct calculation of the critical length we must analyze in detail the mechanisms 
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of growth and shrinkage of a rectangular droplet; these are in general similar to those for 
Glauber dynamics f7|, although the details are different. We will comment below on the 
possibilities of making the following heuristic discussion rigorous. 

Consider first growth. From Figure p] it is clear that a single plus protuberance on one 
of the four sides of the rectangular is not stable; growth proceeds through the formation of 
a double protuberance, which then grows "quickly" (i.e., on the time scale exp 2/3(1 — h)) 
to complete the additional side. The parallel dynamics permits the double protuberance to 
form in one time step, as shown in Figure @(a); the typical time for this process is 

T double ~e 4 ^. (16) 

Alternatively, the protuberance can grow in two consecutive time steps, as shown in Fig- 
ure 0(b); the parallel character of the dynamics enters here as well, since after formation 
of a single protuberance at the first step there must occur, at the second step, both the 
persistence of this protuberance, with probability exp(— 2/3(1 — h)), and the flip of a minus 
spin adjacent to the protuberance, also with probability exp(— 2/3(1 — h)) (see Figure [I]); the 
typical time for this growth process is thus 

r s ^e^xe« (17) 

first step second step 

Clearly r sing i <C r double for (3 large and hence the most efficient growth mechanism is the 
two-step one. 

Again from Figure p] it is clear that the most efficient shrinking mechanism is the usual 
corner erosion, shown in Figure 0(c); the shrinking is performed via a sequence of stable 
configurations. We estimate the time needed for the loss of one of the shorter sides of 
the rectangle, which requires the erosion of Z — 1 sites (after which the remaining single 
protuberance vanishes rapidly). When (3 is large, such a process will typically occur without 
backtracking. The rate at which the entire process occurs is thus estimated as the rate 
for one erosion, exp~ 2/3( - 1+h \ times the probability that I — 2 further erosions occur within 

the lifetime exp 2 ^ 1- ^ of a stable configuration, which is of order 
Thus the shrinking time is estimated as 



exp -2/3(l+h) 2/3(l-/i)l ' 2 



^shrink ~ ° 



2/3(1+/!) x 



e 2j8(l+/0 "I l ~ 2 



,2/3(l-/i) 



18) 



The estimate of the shrinking time r s h r ink can be supported by considering a random 
walk which models what happens on one edge of the droplet. Consider a Markov chain X t , 
t = 0, 1, 2, taking values in the nonnegative integers and with transition probabilities 



P{k,l) 



|e^ c iil — k+1, 



fe-/» ifZ = fc-l, 



otherwise, 
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and 

P(0,1) = ^ C , P(0,0) = l-ie- /3c , 

where k > 1, c > b > 0, and /3 > 0. This chain, with 6 = 2(1 — h) and c = 2(1 + /i), is 
an approximate description of the behavior of the edge of a droplet (see Figure if one 
thinks of X t as representing the number of minus spins on this edge at time t. In order to 
estimate r shrink one should calculate the typical time to see / — 1 minus spins on the edge, 
starting from zero minus spins. Such an estimate, which agrees with fllTf ), is provided by the 
following lemma, whose proof is parallel to that of Lemma 1 of [0. 

Lemma 4.1 For k > 1, define the hitting time r° for the Markov chain X t with X = by 

rl d ^ f {t > 1 : X t = k}. (19) 



Then for any e > 0, 



y f3ck-/3b(k~l)-/3e ^ ^ J3ck-f3b(k-l)+j3e\ P 



To complete the derivation of the critical length for rectangular droplets we compare ([T] 
and ([181) : growth occurs with probability one in the zero temperature limit, that is, / > 
if lim^oo r sing i e /r shrink = 0. This again leads to 



/* 



+ 1. (21) 



We believe that the above argument could be made rigorous along the lines of the cor- 
responding arguments in Section 2 of [0]. The main complicating factor appears to be that, 
because the growth time r sing i e is so large, processes beyond simple corner erosion must be 
accounted for when evaluating the shrinking time. For example, several corners may disap- 
pear in one step; more complicated processes, such as the two-step shrinkage by two sites 
shown in Figure |], are also relevant. These modifications appear to be technical only and 
should not affect the estimate (|T8| ) of r shrink , but we have not carried out a complete analysis. 



5 Monte Carlo results 

The critical length 1^ introduced in the previous section characterizes the behavior of the 
system in the limit (3 — > oo. In this section we define a critical length l% h at finite j3 and 
describe the results of Monte Carlo simulations evaluating 1% h numerically for several values 
of P and h. We find that when (3 is large enough the resulting estimates of 1% h are close the 
theoretical estimate of Z£ given in the previous section. 

Let P/3,h(l) denote the probability that a square droplet of side I grows and covers the 
whole lattice, that is, that in evolving from this initial configuration the system reaches the 
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state +1 before the state —1. Clearly pp,h{l) is a nondecreasing function of I with pp,h(0) = 
and p/3^(L) = 1, so that the differences 

dp, h {l)=Pp, h {l)-pp, h {l-l) ( 22 ) 

form a normalized probability distribution. In the limit (3 —>■ oo (assuming for simplicity 
that 2/ his not an integer), pp^ reduces to a step function, 

= 1; (23) 

and dp^{l) to a unit mass on the critical length At finite temperature, then, we define 
the critical length to be the mean of the distribution dpy. 

£ h = £l<W0- (24) 
i 

From (^) it follows that l* 00 h = l%- We note that this approach in the numerical estimate of 
the critical length is different from that used in ||. 

In the above discussion we have suppressed the dependence of pp^, dp^, and 1% h on the 
lattice size L, since for large (3 and L we expect l% h as defined by (|2~4f) to be essentially 
independent of L. 

We have carried out numerical experiments to estimate the function pp t h(l), and hence 
lp h , for (3 = 0.9, 1.1, 1.3 and h = 0.05, 0.1, 0.2; we varied I over the range of values in which 
Pp,h{l) changes rapidly and made N = 100 runs for each value of I. For each run, we first 
prepared our system in a starting configuration characterized by a single square droplet of 
plus spins of size I, placed in a lattice whose size L was chosen large enough to avoid boundary 
effects. We then followed the evolution of the system and decided, by means of lower and 
upper cutoffs on the total system magnetization, whether the droplet would ultimately grow 
or shrink. We also introduced a cutoff on the total length of each run, chosen as a function 
of (3 so that for most runs the fate of the droplet was determined before the cutoff was 
reached; in the case (3 = 1.3, the highest value of (3 we have considered, this cutoff was 
200,000 iterations. Letting Gp t h(l) denote the number of times that the droplet grew and 
Sp,h{l)t>h.e number of times that it shrank, and assuming that the fraction of the remaining 
iV — G/3 t h — Sp^h cases (in which we did not determine the behavior) in which the droplet 
would have grown if we had waited long enough is the same as for the cases in which the 
behavior was determined, we are led to the estimate 

From the estimated values of pp,h{l) we computed the finite temperature critical length, 
via fl24l) , and the standard deviation of the distribution dp^] the values are recorded in 
Table |T[ The results are in a very good agreement with our theoretical prediction: when 
the temperature is lowered, the numerical measure of the critical length tends to the zero 
temperature theoretical prediction l* h = [2/h\ +1. We did not consider higher values of (3 
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because too long runs would have been needed, but the values we have considered seem to 
be sufficient to see the zero temperature limit behavior. Note that the standard deviation of 
the distribution dp^il) decreases when j3 is increased. This good behavior is clearest in the 
case of small external magnetic field; presumably, higher values of (3 should be considered 
at higher h to approach the limiting behavior. 

We observed that the typical time for growth of the initial square depended strongly on 
j3 and h, but not on I; while the typical shrinking time increase sensibly when I is increased. 



This is qualitatively in agreement with theoretical estimates fllTD and (ITS 



6 Conclusions 

In this paper we have studied the problem of metastable states in probabilistic cellular 
automata, viewing the latter as the simplest instance of models evolving under parallel 
dynamics. All detailed work has been focussed on a particular two dimensional 

model on the square lattice in which the probability of a spin flip at site x depends only on 
the total magnetization of the set of five spins in a cross centered at x (see ([!])). 

We conclude that the general pattern of analysis which has been used for similar models 
evolving under Glauber dynamics applies here as well, since events in which the system 
makes a one-step transition to a configuration significantly different from the current one 
can be neglected in the low-temperature limit. In particular, we argue that the path of 
escape from metastability passes through a critical rectangular droplet. On the other hand, 
the parallel nature of the dynamics does influence the details of the analysis of the escape 
time and, in particular, adds enough complications to make a rigorous analysis more difficult 
than in the Glauber case. 

For the model in question we have shown, through heuristic arguments and Monte Carlo 
simulations, that the critical length of a rectangular droplet is Z£ = \_2/h\ +1. Our theoretical 
prediction is valid only in the limit of zero temperature, but our simulations confirm estimates 
close to the theoretical ones even at finite temperature. 

It is natural to ask whether escape from the metastable state is facilitated or hindered 
by the use of parallel (as opposed to serial) dynamics. It is not clear that this question has 
a universal answer, but as a preliminary approach we may ask what would happen in the 
model of this paper if a serial evolution rule were adopted, so that at each time step one 
spin is chosen at random, with uniform probability, and then updated with probability given 

by (|). 

As remarked in Section |], we used the parallel character of the dynamics only in the 
estimate of r sing i c , so that to estimate of the critical length in the serial case one should 
compare the shrinking time QlSj) with a new growing time r growt h ~ exp[2/3(3— h) +2/3(1 — h)], 
obtained by noting that in the second step of the double protuberance growth the persistence 
probability of the single protuberance need be taken into account. Comparison of these two 
times shows that the critical length in the serial case is given, for h very small, by [3/2/iJ + 1. 
Thus, for these models, the parallel rule leads to a larger critical droplet and a slower exit 
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from the metastable phase. 
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Table 1: Estimates of the critical length l% h , with their standard deviations, obtained from 
Monte Carlo simulations via the procedure described in Section |5]. For (3 = oo we give l* h as 
obtained from fl2~l|) . 
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Figure 1: Probabilities for the flip of the central spin for all possible configurations in the 
5-spin neighborhood. 
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Figure 2: Growth and shrinking mechanisms: (a) double protuberance growth mechanism, 
(b) single protuberance growth mechanism, (c) corner erosion. 
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Figure 3: Approximate description of the behavior of the edge of a droplet. 
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Figure 4: One mechanism occurring on time scale faster than r S i ng i e and hence relevant for 
complete treatment of r s h r ink- 
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